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We study the magnetic properties of graphene edges and graphene/graphane interfaces under the 
influence of electrostatic gates. For this, an effective low-energy theory for the edge states, which is 
derived from the Hubbard model of the honeycomb lattice, is used. We first study the edge state 
model in a mean-field approximation for the Hubbard Hamiltonian and show that it reproduces the 
results of the extended 2D lattice theory. Quantum fluctuations around the mean-field theory of 
the effective one-dimensional model are treated by means of the bosonization technique in order to 
check the stability of the mean-field solution. We find that edge magnetism at graphene/graphane 
interfaces can be switched on and off by means of electrostatic gates. We describe a quantum phase 
transition between an ordinary and a ferromagnetic Luttinger liquid - a realization of itinerant one- 
dimensional ferromagnetism. This mechanism may provide means to experimentally discriminate 
between edge magnetism or disorder as the reason for a transport gap in very clean graphene 
nanoribbons. 

PACS numbers: 73.22.Pr,75.75.-c,71.10.Pm,73.20.-r 



I. INTRODUCTION 

Since graphene can routinely be isolated in the 
laboratoryi^, this monolayer of carbon atoms has re- 
ceived much attention because of its remarkable struc- 
tural and electronic properties. One of the more recent 
graphene riddles is the one about the existence of edge 
magnetism. This phenomenon is based on a simple in- 
tuitive picture: zigzag edges of honeycomb lattices sup- 
port so-called edge states, i.e. exponentially localized 
electronic states at the edges with nearly zero energy. 
The 'flatness' of the energy dispersion of the edge states 
leads to a high local density of states at the Fermi en- 
ergy near the graphene edges. Therefore, these flat bands 
become susceptible to electron-electron interactions; the 
electronic system can lower its total energy, for instance, 
by polarizing the electron spin in the edge states. One 
consequence of edge magnetism is the appearance of a 
transport gap in narrow graphene nanoribbons (GNRs) 
which is inversely proportional to the GNR width. In- 
deed such a transport gap has been measured^. However, 
it cannot be attributed unequivocally to edge magnetism. 
Other mechanisms, like Coulomb blockade in GNRs with 
rough edges*, have been shown to be also compatible 
with the experimental results. 

Though a finite magnetization localized at edges of 
graphene nanoribbons (GNR) has never been directly 
measured, there seems to be a considerable consensus 
on the theory side about the presence of edge magnetism 
in clean GNRs: Hubbard model mean-field theories^—, 
ab-initio calculations^, and even theories that include 
quantum fluctuations^— consistently predict a ground 
state with a finite local magnetic moment at zigzag edges. 
Most of these calculations are numerical and are based 
on two-dimensional calculations on a honeycomb lattice. 
The actual one-dimensional character of edge magnetism 
has been appreciated only by a few authors (see Ref. [H) . 

In this paper, we show that the underlying mecha- 



nism of edge magnetism can be fully understood from 
a one-dimensional point of view, namely by an effective 
model which retains only the edge states while the bulk 
states are dropped. This is a considerable simplification 
but we will show that the deviation from a numerical 
two-dimensional lattice calculation is small. Our effec- 
tive model is based on a simplified edge state model 
which we have used in an investigation of edge states 
at graphene/graphane interfaces^. The crucial feature 
of this model is that it accounts for a finite bandwidth of 
the edge states. The advantages of this effective model 
are remarkable: (a) the mean-field theory is accessible an- 
alytically, (b) the impact of farther neighbor hoppings, 
electrostatic gates, graphene/graphane interfaces, which 
are encoded in the bandwidth of the effective edge state 
dispersion, on the edge magnetism can be investigated 
- also analytically, (c) quantum fluctuations can be in- 
cluded in a large parameter regime within the framework 
of bosonization. 

In addition to this, our theory offers a parameter, the 
bandwidth of the edge states, by which the ferromag- 
netic transition can be driven. We propose a specific 
configuration of electrostatic gates at graphene/graphane 
interfaces which provides direct experimental access to 
this parameter so that the critical regime of this tran- 
sition can directly be investigated experimentally. One 
of the big gains of such an experimental knob by which 
the edge magnetism can be turned on and off is that 
it provides a method to discriminate between different 
sources of transport gaps in clean graphene nanoribbons 
(terminated by graphane): if the transport gap is due to 
edge magnetism, it should disappear as the bandwidth of 
the edge states is increased sufficiently so that the edge 
magnetism is switched off. Disorder induced transport 
gaps, on the other hand, should be largely unaffected by 
this gate because it is designed not to change the carrier 
density in the graphene region. 

This paper is organized as follows. In Section|IT]we dis- 
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cuss the effective edge state modei and the appearance 
of a finite bandwidth on the basis of a single-particle 
picture. The various possible mechanisms which affect 
the edge state bandwidth are subsumed in a single ef- 
fective parameter. Section [TTTI is dedicated to the mean- 
field treatment of the effective electron-electron interac- 
tion in the edge state model. In Section IIVI wc address 
the quantum fluctuations in the edge states by means of 
a bosonization technique. We close with a critical dis- 
cussion of the applicability of our model and the experi- 
mental impact of our findings in Section IVl 



II. GATED EDGE STATE MODEL 

In this section, we cxcmplarily consider an a edge 
in graphene or, equivalently, an a interface between 
graphene and graphanei^. For the /? edgei£, the findings 
are qualitatively similar, however some complications due 
to a possible commensurability with the lattice^ may 
arise. We choose a simplified description of graphene in 
which wc only take into account nearest neighbor hop- 
ping between the it orbitals 
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with t ~ — 3eV. The operator Ci M annihilates an electron 
with spin a at site i. Here, i := (n\,ri2, s), n\,ri2 € 
Z, s = 0, 1 is a collective index which represents the lat- 
tice site at Ri = n^ai + U2&2 + sS. ai, a.2 are the Bra- 
vais lattice vectors and S is the vector which connects 
A (s = 0) and B (s = 1) sublattice sites (see Fig. [T|). 

runs over nearest neighbors on a honeycomb lattice. 
At the edges of the system, the sum must be restricted 
appropriately^. 




FIG. 1: (Color online) Lattice vectors ai,a2 and nearest 
neighbor vector S at an a edge. The coordinate n = 712 iden- 
tifies the position perpendicular to the edge direction. The 
full (green) ellipses indicate the bulk unit cells (n > 1), while 
the dashed (red) circles indicate the truncated unit cell at the 
edge (n = 0). The sublattice index s =A,B is also shown. 

We aim at a description of systems which are lattice- 
translationally invariant along the a! direction and it is 
convenient to transform this direction to k-spacc. There- 
fore, we work with the electron operators d n k s a — 



N x 2 J2 ni e lknic ni,n,s,a, henceforth. N x is the number 
of unit cells in the ai direction. It is well known22 that 
in this model, a zero-energy state exists which is expo- 
nentially localized at the a edge 
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where = —(In \uk\) 1 is the localization length, A4 
: 1 — |it fc | 2 )2 is a normalization constant, and Uk = l + e lk . 
<f> is some unimportant phase. Note that this a edge state 
exists for k € 4f ] , the rest of the Brillouin zone being 
the domain of the (3 edge stated. 

We have shown in Ref. [IH that it is possible to incor- 
porate some details of the edge states, like a finite band- 
width or spin-orbit interaction, in this simplified model 
on an effective level. The terms in the Hamiltonian which 
are responsible for these additional edge state properties 
are usually much smaller than Hq. Thus we may treat 
these small terms within perturbation theory by project- 
ing them onto the subspace which is spanned by the edge 
states \tpk,(r}- 

The dominant effect of a certain class of Hamiltoni- 
ans describing, e.g. graphene/graphane interfaces^ or 
next-nearest neighbor hoppings^, is to create a non-zero 
bandwidth of the edge states. We model this class of 
bandwidth-generating edge properties as an effective lo- 
cal gate, described by the Hamiltonian 
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where t e is the amount by which the on-site energy of 
the outermost carbon atoms of the a edge are changed. 
Note that t e is an effective parameter which comprises, 
apart from the potential due to a real electrostatic gate, 
also many other details of the edge. As long as |t e | <C \t\, 
we may resort to first order perturbation theory in Hq 
and find for the energy of the a edge state \4>ko) 



eo(fc) ~ — £ e (2cos(fc— 7r) — 1) 
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Usually an edge gate is not atomically sharp. However, 
as long as the gate is localized at the edge, the qualitative 
appearance of the dispersion is insensitive to the 'leakage' 
of the gate potential into the graphene bulk. 

A mechanism based on graphene/graphane interfaces 
which admits an especially well tunable dispersion of the 
form (j4|) is discussed in Appendix [A] The tunability of 
the bandwidth t e turns out to be crucial for changing the 
magnetic state of the edge. t e can be on the order of eV. 

At this point, we would like to give some arguments 
why it is sufficient for the analysis of edge magnetism to 
keep only the edge states while the bulk states are com- 
pletely removed from the considerations. These heuristic 
arguments arc supplemented by a direct comparison be- 
tween the effective model and a full lattice calculation in 
mean- field approximation (see Appendix IXT|) . 
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FIG. 2: (color online) Non-interacting band structure of a 
narrow a/3-ribbon (20nm wide) with graphane terminations. 
The calculation is based on a tight-binding model which takes 
into account the ir- and cr-orbitals of the carbon sitesi^. The 
dashed line indicates a typical Fermi energy (see text). 

The most straightforward argument is based on the 
band structure of a/3 ribbons^. At reasonable fillings (as 
shown in Fig. [2]) only the edge states cross the Fermi en- 
ergy while the bulk states are energetically remote. Since 
the a and /? edge states are localized at their respective 
edges, their mutual spatial overlap is exponentially small 
as long as k does not get too close to one of the Dirac 
points K,K'. From the combination of the energy argu- 
ment and the localization argument follows that the a 
edge state can be examined independently of all other 
states. The same argument holds also for the /3 edge 
state. 

How is the argument to be changed if, instead of aft 
ribbons, we consider the usual aa ribbons? Then, in- 
stead of having a /3 edge state attached to the a edge 
state at K and K', the edge state merges into a bulk 
state which is not exponentially localized at the other 
edge. However, the overlap between the edge state and 
this bulk state is still small (though not exponentially 
small) because the bulk state wave function is essentially 
proportional to W~? s'm(iry/W), with W the width of 
the ribbon, and thus is small where the edge state wave 
function is large. 

Finally, one may ask what happens if we drop the re- 
striction of narrow ribbons. In this case, the energy ar- 
gument fails because it relied on the presence of a finite- 
size gap. Now, we resort to a rough scaling analysis of a 
density-density interaction like the Hubbard interaction, 
which is used below, or any screened interaction. In such 
a type of Hamiltonian, the spatial overlap of the wave 
function weights is important. The spatial density over- 
lap of an edge state with localization length £ with the 
nth bulk state is roughly 

o{n) ~^^dysm 2 (nny/W) ~ (5) 
For increasing W, there are ~ W bulk states at the Fermi 



energy so that the sum of all density overlaps of the edge 
state and the bulk states is Y^Z^i ~ i- e - tne 
total overlap does not grow with W. Thus, we expect 
that the bulk states do not become important as W — > 
oo because the spatial overlap between a typical edge 
state and the energetically relevant bulk states decreases 
sufficiently fast with W. 



III. HUBBARD INTERACTION IN THE 
PROJECTED MODEL 

We want to study the electron-electron interaction 
effects in the edge states in the Hubbard approxima- 
tion Hjj = U ^\ PtfPii with the local density operator 
Pi a = c[ a Ci a . In terms of the d n ^k,s,a operators we have 

Hu = UN~ 1 Pn,q,8,tPn,-q,s,U (6) 

n,q,s 

With p n . q .s,a = T,k d i,k+ q ,s^ d n.k,s,a- 



A. Projection onto the effective edge state model 

Although the interaction energy scale U can be quite 
large, we start with considering Hjj as a perturbation 
to Hq . We show in Appendix |D] by a comparison to a 
less restricted numerical calculation that the following 
perturbative calculations capture the essential physics. 

Since we aim at an effective description in terms of 
the edge states, disregarding the bulk states, we need to 
account for the background electron density, generated 
by a completely filled valence band. In the bulk, charge 
neutrality is reached if the complete valence band is filled 
and the conduction band is empty, i.e. one half of the 
7r band is filled. At half-filling there is on average one 
electron per site in the bulk system. In the bulk one can 
easily convince oneself that, because of the A B sublat- 
tice symmetry^, there really is exactly one electron per 
site - not only on average. At an edge, the sublatticc 
symmetry is broken so that the site occupation may dif- 
fer. Here, the edge states come into play. Half-filling is 
reached by occupying the valence band and one half of 
the edge states. The answer to the question, which half of 
the edge states have to be filled, depends on the energetic 
details of the edge states which are due to, e.g., spin-orbit 
interaction, dispersion and electron-electron interaction. 
If, for instance, we assume a non-magnetic configuration 
such as the one considered in Ref. [l9| in which the spin- 
up states for k G it] and the spin-down states for 
k G [n, 4^] are occupied, it turns out that also at the 
edges each site carries exactly one electron charge - one 
half with up-spin and one half with down-spin. The con- 
tribution to the total electron density per spin at unit 
cell n and sublattice site s, coming from edge states that 
are occupied in this way is (note the invariance of the 
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edge state wave function under k — > 2n — k) 



* dfc 

JL 2^ 



Nt[2{l + vo S k)] n =:-p {n,s), (7) 



where po(n,s) is the electron density per spin due to a 
fully filled edge state band. The background electron 
density Pb(ti, s,a) (relative to half-filling) of a filled va- 
lence band and an empty edge state band is therefore 
given by 



p B {n,s,<r) = --p (n,s). 



(8) 



A more rigorous derivation of Eq. ([5]) can be found in 
Appendix [CJ 

In the remainder of this work, we will drop the sub- 
lattice index s, setting it to B. This is because the A 
sublattice sites are always half-filled in the perturbative 
treatment and thus do not contribute to the following 
considerations. 



B. Polarized vs. unpolarized state 

In order to gain a rough overview over the basic princi- 
ples from which the magnetic properties of the graphene 
edges derive, we start with comparing the energy of 
a ground state configuration with fully polarized edge 
states \P) = rj fc 4,t to thc 

energy of an unpolarized 
configuration \U) = H^Hsjl^i^ e\ a \9) . |0) repre- 
sents thc completely filled valence band, which plays thc 
role of the vacuum in the effective model. The opera- 
tor el creates one edge state with crystal momentum k 
along the edge and spin a. 

The total kinetic energies (per unit length) of these 
states, relative to the kinetic energy of a completely filled 
valence band, are 



! 3 (9) 



t e E v := (U\H , G \U) = -2t e 
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where the numerical factors Ep and Eu have been de- 
fined for later convenience. The Hamiltonian of the ki- 
netic energy Hq q = Hq + Hq consists of the usual hop- 
ping Hamiltonian of graphene Hq and of the gate Hamil- 
tonian Hq, which is responsible for thc finite bandwidth 
of the edge states. 

For the calculation of the Hubbard energy, we treat 
all electron densities relative to half-filling. In the fully 
polarized state \P) the spin-dependent electron density 
is given by 



pp{n,a) = cr-po(n). 



(11) 



In the unpolarized state | U) , the electron density is spin- 
independent 



Pu(n,a) 
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From the electron densities, the interaction energies are 
given by U^n^o Px(n,t)px(n,i), where X = P,U. 
This leads to 



UW P := -U 



~ dfcdfc' 
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for the interaction energy of the polarized state and 
^dfcdfc' s k s k ,{N k N k ,) 2 
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for the interaction of the unpolarized state. 

Obviously, for fiat bands (t e = 0), thc polarized state 
\P) is thc ground state. If the bandwidth t e is increased 
to positive values, however, the total energy of the polar- 
ized state Ept e + WpU eventually becomes larger than 
the total energy of the unpolarized state. This happens 
at 



W P - W v 
Ep — Ejt 



0.17, 



(15) 



the critical bandwidth/interaction ratio for positive t e at 
which thc system becomes unpolarized. For negative t ei 
the polarized state becomes unstable with respect to an- 
other unpolarized state that has the inverse edge state 
occupation of \U). This instability occurs at the sym- 
metric position t e /U ~ —0.17. 



C. Solution of the mean- field equations 

The argumentation in the previous paragraph is of 
course very superficial since we have only compared the 
total energies of completely polarized and completely un- 
polarized states, neglecting thc possibility of partial po- 
larization. Thus, we now formulate a closed set of mean- 
field equations for thc effective edge state model with 
Hubbard interaction and solve them. 

At zero temperature, the mean-field energy of an edge 
state \ipka)> with all other edge states with mean-field 
energies smaller than ep occupied, is given by 
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(16) 
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where we have defined the effective interaction vertex 



r(k,k',q)= k+q k k q k 



1 



l k + q U kU k >- q U k' 



(17) 



r(fc, fc', 0) results from the summation of the probability 
densities of two edge states, k and fc', over the sites in a 2 
direction. For later convenience, T has been introduced 
in a more general form than actually needed here, namely 
for non-zero q. It is important to note the invariance of 
r(fc, fc', 0) under fc — s- 2n — fc and fc' — s- 2tt — fc', as well as 
under fc o fc'. 
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FIG. 3: (Color online) Mean-field result for the edge mag- 
netization Me as a function of the bandwidth t e of the edge 
state. The topmost, thick (red) line represents the result for 
half-filling. The other lines show the results for increasing 
the filling from 1/2 (topmost line, red) to 0.85 (bottom line, 
blue). The dashed vertical lines indicate the critical band- 
width obtained from the energy comparison of polarized and 
unpolarized states at half-filling. 



The Fermi energy is fixed by a constant effective elec- 
tron density n. n = means that the edge states are 
completely unoccupied and n = 1 means that all edge 
states are occupied. 
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(18) 



Eqs. (|T6|) and (jT8|) arc a complete set of mean-field equa- 
tions which can easily be solved. It is convenient to char- 
acterize the solution by the edge magnetization which we 
introduce as 



M, 



= 5> P (0,(7) 



(19) 



where p(n,a) = J f^6(e F - e CT (fc)) |(0| d n ,k,B,o \ipk,a)\ 2 is 
the spin-dependent density at site n. 

Fig. |3] shows the solutions of the mean-field equations 
for different fillings n. One observes a second order quan- 
tum phase transition at a critical \t e /U\ CT it_, below which 
the edge becomes spontaneously polarized. The polar- 
ization saturates as \t e /U\ — > 0. For half-filling, the satu- 
rated regime corresponds to the usual edge magnetism^— 
which is characterized by a complete polarization of the 



edge state spin. If h > |, also the minority spin edge 
state becomes partially occupied until, for n = 1, both 
spin species are completely filled and no spin-polarization 
exists any longer. 

Note that, in addition to the phase transition at which 
the polarization starts to grow from zero, there is also 
another transition at which the polarization starts to de- 
viate from the saturation value. While the first transition 
is of a Stoner type, the latter is not. As an illustra- 
tion of the mean-field solution we provide three movies^ 
showing the single-particle energy e cr (fc) (see Eq. (|TB)) ) 
for both spins together with the edge magnetization for 
three different filling factors h = \ , | , § . It can be seen 
that the situations n = j , | are connected by particle- 
hole symmetry. Furthermore, we see that, at the Stoner 
transition, the Fermi points of the two spin species start 
to become more and more split . At the second transition, 
which is not of Stoner type, the number of Fermi points 
changes from two to four. 

In the following, we distinguish between the saturated 
regime and the regime close to [t e /U] clit where the spin- 
polarization is small. The latter is called the regime of 
weak fcrromagnctism. This regime will be especially im- 
portant for the analysis of quantum fluctuations. Note 
that the stability of weak ferromagnctism on the mean- 
field level is a consequence of the special momentum- 
dependent form T(k,k',q) of the effective edge state 
model. For a truly one-dimensional Hubbard model, the 
interaction vertex would be constant in fc, so that a so- 
lution of the self-consistency equations does not yield 
weak ferromagnctism. We note that the mechanism we 
describe here is different from the one Bartosch et al. 
use^i for the stabilization of weak ferromagnetism in one- 
dimensional metals. 



D. Weak ferromagnetic regime 

The general mean-field equations have been solved nu- 
merically in the previous subsection. Close to the critical 
point, however, one can obtain an approximate analytical 
solution with 



k 

u 



(20) 



as a small parameter. This will be useful below where 
quantum fluctuations around the mean-field solution in 
the weak ferromagnetic regime (but not too close to the 
transition) shall be analyzed by means of the bosoniza- 
tion technique. 

Weak ferromagnetism is characterized by a small im- 
balance in the spin population, which is quantified by 
spin-dependent Fermi momenta fc^ro- = tt + r(kp + trAfc) 
(see Fig. @|. Here, r = R(L) stands for right- moving 
(left-moving) parts of the edge state dispersion at the 
Fermi energy. Henceforth, R = +1 and L = — 1 when 
used in equations. Likewise, a = ±1 stands for up- 
and down-spin, respectively. We use the symbol kp for 
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FIG. 4: Fermi levels in the non-magnetic ground state and 
in the magnetic ground state where the Fermi momenta kFra 
are split by 2Ak in k-direction. 
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FIG. 5: Critical bandwidth [t e /U] C Ht. as a function of the 
edge state filling. stands for a completely empty edge band 
(fc.Fr = 7T in the non-magnetic configuration) and 1 for a com- 
pletely filled edge band (&Fr — it + rir/3 in the non-magnetic 
configuration). 



the distance of the two Fermi points from ix in the non- 
magnetic configuration, where Afc = 0. A fixed k F cor- 
responds to a fixed edge state filling n, in both, the mag- 
netic and the non-magnetic ground state. 

For non-zero Afc, the kinetic energy and interaction 
energy are 



A A 



dfc 

— sign(fc)e(7r + k F + k) (21) 
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r(fc t ,fc i; 0), (22) 



where we have defined the single particle kinetic energy, 
corrected by the background charge density of the filled 
valence band, e(fc) = erj(fc) + Uns(k), with 



n B (k) 

We approximate for small Afc 
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(23) 

(24) 
(25) 



where we have dropped terms of order Afc 3 in the kinetic 
energy and terms of order Afc 5 in the interaction energy. 
The 0(Afc 3 ) terms in the kinetic energy could have been 
taken into account but it turns out that they only lead to 
an inessential quantitative rcnormalization of Afc. This 
approximation is equivalent to linearizing the single par- 
ticle dispersion e(fc) ~ Vp(k — k Fra ) around the Fermi 
points. 

In the ferromagnetic regime, the total energy assumes 
a minimum for 



Afc~± 



6a 
— -; 

73 



for t > 0. 



(26) 



If we had taken the third order terms in Eq. (f24|) into 
account + ^), the right hand side of Eq. ([2"u]l 



would have been multiplied by a factor of ~T.02. Thus, 
it is reasonable to neglect the curvature of e(fc). 

The parameters a, /?, 71 , 73 can be calculated from Eqs. 
pi) and (|22|) . This is especially convenient for a 'magic' 
filling, corresponding to kp = arccos|, which is slightly 
lower than half-filling. Throughout the remainder of this 
paper, all analytic discussions are based on this filling. 
The reason for this is simply that the formulas are only 
about one third as long as for half-filling, for instance. 
We find 
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(27) 
(28) 

0.209 (29) 
~ 1.96 (30) 



and thus a critical bandwidth 



u 

= 2 



71-/3 



V3(175 - 68a/5) + 200tt - 600 arcsin | 

V ; §- ~ 0.244. 



75tt 



(31) 



We note that these parameters can be calculated ana- 
lytically also for general fillings. However, the results are 
quite cumbersome. Wc therefore Fig. [5]shows the critical 
bandwidth as a function of the edge band filling. 

Note that we describe essentially a Stoncr mechanism 
here. This is why we will call the critical point r = the 
Stoncr point, henceforth. 
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IV. QUANTUM FLUCTUATIONS 

In the previous section, we have analyzed the mean- 
field theory of the effective edge state model. Now, we 
proceed by analyzing the stability of the mean-field so- 
lution with respect to quantum fluctuations. Since the 
effective model is one-dimensional, quantum fluctuations 
may be treated most easily by means of the bosonization 
technique^. We have shown in the preceding section that 
the linearization of the single-particle dispersion around 
the Fermi points is a good approximation. This is an es- 
sential prerequisite for the applicability of bosonization. 

With the annihilation operator e^ a of an edge state 
with momentum fc along the edge and spin er, the full 
Hamiltonian of the effective edge state model can be writ- 
ten as 



= Y J Mk) + Un B {k)]e\ a e ka 

k.a 



— r (M^)4+?t efc t4'- 9 i e fc'i> ( 32 ) 



N. 



k,k' ,q 



where the primed sum means that the summation is re- 
stricted such that ^ < fc, k', k + q, k' — q < 4p For the 
derivation of the interaction term, see Appendix [B] The 
background term ns(k) describes the additional disper- 
sion coming from the electron density of the filled valence 
band. We introduce a normal ordering 



A: = A — (<j) \A\(j>o) , 



(33) 



where \<fio) denotes the Slater-determinant of the ground 
state of the mean-field approximation to H, as discussed 
in the preceding section. With the generalized density 



n a (k) = J2T(k,k',0) (4> \e 



the mean-field part of H becomes 



.efe'er 



(34) 



H a = 



k.a 



eo(fc) + U(n B (k) + n_ CT (fc))] e\ a e ka (35) 



and the term describing the quantum fluctuations around 
the solution of reads 



U 



zr r (k,k',q) : 4+ q t efc t :: 4'-?i e fc'l : • ( 36 ) 

k : k' ,q 



Note that the mean- field energy e CT (fc), defined in Eq. 
([33]) is consistent with Eq. ([TB"]) . 

For the bosonization we need to distinguish between 
the non-magnetic mean-field phase, which will turn out 
to behave as an ordinary spinful Luttinger liquid, and the 
magnetic phase, the properties of which are somewhat 
more intriguing. Especially the boundary of these phases 
will turn out to be complicated so that, in this work, we 



restrict the discussion to values of t e /U ', sufficiently far 
from the critical point. However, t e /U is still required to 
be larger than the value at which the spin-polarization 
of the edge-states saturates. 



A. Non- magnetic phase 

In the non-magnetic phase, i.e. for t e /U > [t e /U] CTit , 
both spin species are equally occupied so that e a (k) is 
spin-independent. The Fermi velocity is given by 



2t e sin(kp) 



-U 



^ dfc 
! 2^ 



e{k F - \k\) - 



dfc, 



-r(7T + fcF,7T + fc,0). 



(37) 



An essential simplification needed in order to express H 
by means of bosonic degrees of freedom is the lineariza- 
tion of the edge state dispersion e a (k) around cf at the 
left-moving (r = L) and right-moving (?* = R) Fermi 
point, i.e. 

e a (k) ~ v F r{k — k Fra ), for \k — fepro-l ^ hp- (38) 





FIG. 6: Mean-field configuration of the non-magnetic regime. 
Both spin species are equally occupied (indicated by the bold 
black lines). 

At the magic filling, to which we want to restrict the 
discussion, the Fermi velocity is 



with 



vf = t e 



15 



US 



_3_ 

5tt 



0.046 



(39) 



(40) 



Because in the Hubbard model only densities of dif- 
ferent spin projections interact with each other one finds 
tha1>2£ gq = 0, Vi. For gi±, the g-ology for the Hub- 
bard interaction in the edge state model (see Appendix 
IB"]) gives 



9l± = 92± = 9i± = Ur(ir + k F , 7r + k F , 0) 



3J7 
"5"' 



(41) 
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Following the standard procedure of Abelian bosoniza- 
tion (see, e.g., Ref. [l8l ). one finds 



H = H° + H° + Hi 



(42) 



where the free Hamiltonians H® and H® for the charge 
and spin sector, respectively, are 



H° v = — / dx 

v 2tt 



u v K v {d x d v (x)f + -^(d x 4>„{x)) 



A" 



where v — c, s with 



(43) 



u c 3U , . 

u c K c = v F and — = 1 + (44) 

K c birvF 



for the charge sector and 



IT A Us 1 iU 

u s K s — v p and = 1 



(45) 



for the spin sector. The bosonic fields 4> v (x) and Q v {x) 
satisfy the commutation rules 

\4> v {x\d x Q v ,{x')\ =md v „>8{x-x'). (46) 

The backscattering process g\± leads to a sine-Gordon 
term in the spin sector 



dx cos 



(47) 



where -q ~ kp 1 is an ultraviolet cutoff^. 

Due to spin-charge separation, the spin and charge sec- 
tors may be treated separately. We want to focus on the 
spin sector here. The spin velocity u s becomes singular as 
5 ^ — > 1. This is exactly the Stoner point of the mean- 
field theory. Note that the Stoner point is also reflected 
in the spin susceptibility of the free bosonic theory, i.e. 
the theory in which the backscattering term Hi is disre- 
garded, xo = 2^T ^ \5-kvf-3U\ ■ X° has a singularity at 
t = 0. 

It is well knowr>i& that a perturbative treatment of 
Hi leads to a renormalization of K s and the amplitude 
of the backscattering process <?ij_. The g\ j_ process is 
marginally irrelevant for SU(2) invariant systems like the 
edge state in the non-magnetic regime. This means that 
K s — > K* = 1 and g±± —> 0, as the ultraviolet part of 
the Brillouin zone is integrated out successively. In the 
rcnormalizcd theory, the spin susceptibility becomes 



Xo 



K* 

27TW, 



1 



oc 



y/\5T7V F - SU\ 



(48) 



However, the renormalization group for the backscatter- 
ing process is perturbative in gi±/u s and cannot be used 
close to the Stoner point, where u s — > 0. This is why 
we exclude the critical region from our argumentation in 
this work. 



B. Weak ferromagnetic regime 

In the previous subsection, we found that the energy 
which has to be paid for fluctuations of the field 4> s (x) 
is proportional to (1 — 3U /5-KVF){d x <p a {x)) 2 . The sign 
of this term becomes negative beyond the Stoner point. 
This means that the system can lower its energy by de- 
veloping spin fluctuations and becomes unstable against 
a new ground state with a spontaneously broken sym- 
metry. However, the bosonic theory based on Eqs. (|43|) 
- (|47jl is not able to actually predict the proper ground 
state. Therefore, we need to go back to the level of the 
mean-field theory. 

We have shown in the previous section that the mean- 
field theory becomes unstable with respect to a spin- 
polarized ground state at the Stoner point. Thus, it is 
reasonable for t e /U < [te/U] crit to bosonize the quantum 
fluctuations around the ferromagnetic mean-field theory, 
rather than the non-magnetic. The question to be an- 
swered is then, if this polarized mean-field theory is sta- 
ble with respect to quantum fluctuations. 




AkAlc 



AkAk 



FIG. 7: Mean-field configuration in the weak ferromagnetic 
regime. The fat lines represent the occupied states and the 
thin lines represent the unoccupied states. The red arrows 
indicate the backscattering process, which is not momentum- 
conserving here. 

We restrict the discussion to the regime of weak fer- 
romagnctism, i.e. we exclude the saturated regime. It 
should be noted that the different spin species may have 
slightly different Fermi velocities. This leads to a term 
(vpf — vpi)(6' s 8' c + <j)' s (j)' c ) in the bosonized Hamiltonian. 
This term mixes the spin sector and the charge sector. 
However, vf^- — Vfi is of the order of Afc which is small 
in the weak ferromagnetic regime. Thus, we neglect this 
difference in the remainder of this work, keeping only the 
mean- value vf = o( VF t + v Fl)- We begin with the am- 
plitudes of the forward scattering processes which lead 
to the quadratic part of the bosonic Hamiltonian. In the 
non-magnetic mean-field ground state, these two ampli- 
tudes are g-2± = go. = UT{it + fei?,7r + fci?,0). In the 
ferromagnetic mean-field ground state, characterized by 
a finite Afc, however, the amplitudes are 

92± = 9i± =■ UT Q (Ak) 

= Ur(TT + k F + Ak,TT + k F -Ak,0) (49) 

To (Afc) has a maximum at Afc = 0, which means that 
the interactions between the electrons at the Fermi level 
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become weaker as the polarization ~ Afc becomes larger. 
On the mean-field level this leads to the balance be- 
tween the kinetic and the interaction energy which allows 
the existence of the weak ferromagnetic regime. In the 
framework of bosonization, the polarization dependence 
of the interaction has the important function of restrict- 
ing u a /K s to positive values: for small Afc we find (for 
the magic filling) 

r (Afc)^^Afc 2 (50) 

so that for t > 

Us/vf ( 136a 5tt fb\ 

— ^{i^--rz) T ^ +i - bT - (51) 

Thus, the free bosonic theory in the ferromagnetic regime 
is meaningful if the proper mean-field theory is used. 




0.15 0.2 0.25 0.3 0.35 

tj U 



FIG. 8: The prefactor u s /K s of the term (d x (j>s) 2 in the 
bosonized Hamiltonian near the Stoner point, as a function 
of t e /U for varying fillings near half-filling. The red curve 
stands for n — 0.4 and the blue line stands for n = 0.6. 

Fig. [8] shows that the previous considerations are not 
special to the magic filling: for all fillings of the edge 
state, u s /K s is non- negative on both sides of the Stoner 
point, if the proper mean- field theory is used. 

Finally, we consider the backscattering term in the 
ferromagnetic regime. From Fig. [7J one can see that 
backscattering is not momentum-conserving for Afc ^ 0. 
However, the momentum mismatch is 4Afc so that close 
to the Stoner point, the backscattering term is at least 
approximately momentum-conserving. A careful deriva- 
tion of the backscattering term gives 

H\ = t^T2 / dxcos (y%M x ) + 4Afca;) . (52) 

For Afc sufficiently large, H\ averages to zero because 
of the fast oscillation with x in the cosine; this is just 
another way of saying that backscattering is momentum 
non-conserving and thus not allowed. 

For small Afc the situation is dramatically different, 
though. Since the prefactor of the free Hamiltonian is 
very small near the Stoner point, H\ is dominant. Thus, 



the field 4> s (x) becomes locked to — 4Afcx. This means 
that the finite mean-field magnetization around which 
we have expanded the quantum fluctuations is brought 
back to zero by the backscattering term, and this violates 
self-consistency. Just as in the non-magnetic regime, our 
treatment becomes invalid near the Stoner point. 

Deeper in the ferromagnetic regime, for r sufficiently 
large, the free Hamiltonian becomes stronger again so 
that the zero-magnetization state becomes energetically 
unfavorable. There, Hi is not strong enough to desta- 
bilize the ferromagnetic mean-field theory and the fer- 
romagnetic ground state is self-consistent. This is 
quite what one would expect from quantum fluctuations, 
namely that they reduce the tendency towards a broken 
symmetry ground state. The critical point, at which the 
system polarizes spontaneously, is shifted to higher inter- 
action strengths. 



V. DISCUSSION 

In the preceding sections we have demonstrated that 
edge magnetism at graphene/graphane interfaces can be 
understood on the basis of a one-dimensional effective 
model for the edge states. This effective model comprises 
a non-magnetic regime for a sufficiently large edge state 
bandwidth, a weak ferromagnetic regime for intermedi- 
ate bandwidths and a saturated regime for small band- 
widths. The saturated regime corresponds to the usual 
edge magnetism as it has been studied in Refs. [5T4l2l 
for instance. The underlying mechanism is based on the 
one-dimensional version of the well known Stoner insta- 
bility. The essential difference to previous works is that 
we utilize the bandwidth t e of the edge states in order to 
obtain a well controlled theory. At graphene/graphane 
interfaces, t e is experimentally accessible by means of 
electrostatic gates, so that the edge magnetism can be 
turned on and off dynamically in an experiment. 

Such experimental control over the magnetic state of 
graphene edges might prove useful for distinguishing be- 
tween edge magnetism and other sources like disorder as 
the underlying mechanism for transport gaps in GNRs. 

The advantage of this control parameter for the the- 
ory becomes evident by comparing this work to Ref. Il2l . 
where the linear kinetic energy term is absent. There 
the exponential overlap between adjacent edge states is 
used in order to define the (non-interacting) kinetic en- 
ergy term around which the quantum fluctuations can be 
treated. It is obviously hard to study the effect of quan- 
tum fluctuations directly in the fully polarized regime 
where the band of one spin species is either completely 
occupied or completely unoccupied, because it is not pos- 
sible to bosonize this theory in the usual way. In our 
effective edge state model, instead, we can tune from the 
usual Luttinger liquid to the ferromagnetic region in a 
well controlled manner. In the intermediate regime of 
weak ferromagnetism (not too close to a certain critical 
point, though), bosonization works well and we are able 
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to study the effect of quantum fluctuations in a ferro- 
magnetic Luttinger liquid. 

It should be noted that ferromagnetic Luttinger liq- 
uids have been studied by Bartosch et al£^. There, it has 
been shown that the ferromagnetism can be stabilized by 
a large positive third derivative of the dispersion. Edge 
states in graphenc or at graphene/graphane interfaces 
do not fall into this category, however. In graphene edge 
states, the ferromagnetism is rather stabilized by means 
of effective velocity-dependent electron-electron interac- 
tions: the larger the distance between the Fermi mo- 
menta of the different spin species, the smaller is the 
effective interaction. 

The velocity-dependence of the interaction comes from 
the momentum-dependence of the localization length of 
the edge state wave function. Therefore, the same mech- 
anism that gives rise to weak ferromagnetism at graphenc 
edges, might be found also in other systems with similar 
localization properties of edge states, e.g. in topological 
insulators where the edge states also become more and 
more delocalizcd as they merge into the bulk. 

Since we are investigating magnetism in one dimension, 
a comment about the applicability of the Lieb-Mattisi^ 
theorem, which forbids ferromagnetic order in one di- 
mension, is mandatory: this theorem is not applicable 
to edge magnetism because the assumptions of Lieb and 
Mattis exclude velocity-dependent interaction potentials. 
Indeed this velocity-dependence is essential for the sta- 
bilization of the weak ferromagnetism. Nevertheless, the 
statement of stability is only a statement of local stability. 
Within the line of argument of the present work, we can- 
not exclude the presence of another phase with lower en- 
ergy than the ferromagnetic phase. More effort is needed 
for a final answer to this question. However, our effec- 
tive model provides a convenient framework for further 
investigations. For instance, it allows numerical quan- 
tum many-body simulations with fewer effort than for 
the full two-dimensional lattice model because all unim- 
portant degrees of freedom have been eliminated already, 
while the important properties of edges of honeycomb 
lattices have been condensed into a 'small' effective one- 
dimensional model. 

However, there is a further indication which helps 
to gain more confidence in the existence of the weak 
ferromagnetic regime: in the limit t e — ¥ 0, the two- 
dimensional honeycomb lattice model becomes particle- 
hole symmetric so that Liebs theorem^ applies and pre- 
dicts a high-spin ground state at a zigzag edge. The 
particle-hole symmetric case corresponds to the com- 
pletely saturated regime in the present work. This 
means that our prediction of weak ferromagnetism at 
graphene/graphane interfaces is consistent with the well 



accepted limits (a) edge magnetism in graphene for t e — > 
and (b) Luttinger liquid behavior for U <C \t e \. 

Another important issue, connected to the low dimen- 
sionality, is the question of the stability of edge mag- 
netism w.r.t. finite temperatures. Of course, the usual 
entropy argument, which prohibits a spin-polarization in 
a one-dimensional system in the thermodynamic limit, is 
applicable here. However, graphene structures are usu- 
ally of mesoscopic size, which means that spin-waves ex- 
hibit a finite-size energy gap. Thus, the ground state, 
which we have studied in this work, should be observable 
for sufficiently small mesoscopic structures at sufficiently 
low temperatures. 
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Appendix A: Graphane gate 

An atomically sharp gate as described above is exper- 
imentally not feasible. Even if the termination of such a 
metallic gate could be controlled on an atomic scale, the 
potential would leak into the bulk graphene region be- 
cause of the weak screening in graphene, thus changing 
the electron density in the graphene region. 

To circumvent this, we propose to use gates of differ- 
ent voltage on different sides of the plane. This gate 
lifts (lowers) the on-site energy of the hydrogen atoms 
on the +z (—z) side of the z = plane. Since the 
graphene atoms are located at z = 0, the on-site energy 
is unchanged in the graphene region. We use the effec- 
tive model of edge states at graphene/graphane interfaces 
from Ref. [l6| which is defined by the Hamiltonian 

/ o m k \ 

H eS = [Wk f , (Al) 
V t' e H ) 

where t ~ —3 eV is the usual hopping between carbon 7r 
orbitals, t' ~ — 6 cV is the hopping between the carbon 
7r orbital and the hydrogen Is orbital in the graphane 
region, and eu — —0.4 eV is the hydrogen Is orbital 
energy without any gates. Since eu is smaller than all 
other energy scales, we treat it in perturbation theory. 
Setting eu = 0, we obtain a zero energy eigenstatc 
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m(k))+M£-\H). (A2) 
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From Eq. (|A2[) . we see that the wave function of the 
edge state at a graphene/graphane interface lives pre- 
dominantly at the graphene region sites where the usual 
graphene edge state lives (|^o(&))) an d at the hydrogen 
atom in the first graphane row. Obviously then, a non- 
zero on-site energy at the hydrogen atom leads in first 
order perturbation theory to an energy shift of the order 
of this on-site energy. The hydrogen on-site energy is 
composed of the intrinsic chemical potential of hydrogen 
en and the gate-induced on-site potential A. Thus, we 
find 

e a (fc) = (e H + A)^(2cos(fc-^)-l). (A3) 

The distance between the hydrogen atom planes 2.8A 
and the distance between the carbon planes is 0.6A. 
Thus, if the hydrogen atom is at the potential A then the 
carbon atom attached to this hydrogen atom is roughly 
at the potential 0.2 A. This gives an additional positive 
contribution to the edge state bandwidth (ejj + A)t 2 /t' 2 , 
since the edge state wave function has also a weight of 
order ch /t on this carbon atom. 

The largest electric fields that can be reached are of 
order 10 7 — 10 8 V/cm. Thus, the total hydrogen on- 
site potential can be of the order of one eV, leading to 
an in-situ tunable bandwidth range of a several hundred 
meV in addition to the edge state bandwidth contribu- 
tions coming from, e.g., farther neighbor hoppings. Also, 
the bandwidth can be tuned by substituting the hydro- 
gen atoms in the graphane region with other elements or 
molecules with different orbital energies. 

Appendix B: Effective Hubbard interaction 

The edge state operator reads 

oo 

e fe , CT =Nk^[-u* k } n d n ^ B ,<T (Bl) 

n=0 

(we drop the sublatticc index henceforth because the edge 
state lives on the B sublatticc only) and its commutation 
relation with the d-opcrator is 

=^4H4]"<W<W- (B2) 

Now we project the Hubbard Hamiltonian (restricted to 
B sites) 

h u = un- 1 ]T 4+^,A»,t4'- 9 ,n,A>'a (B3) 

k,k' ,q,n 

to the Fock space spanned by the edge states (e^ is a 
short form for eki.a) 

f 

Hu= (lt,2i|ffa|3t,4 4>el t 4 4 e4+e 3t , 

(B4) 



where the two-fermion states are [ltr, 2a') — e[ a e' 2a , |0). 
The primed sum is restricted to values of ^ < ki < 
We calculate the matrix element 

Jf Y (°l e a4- e lt4+9,n,t d *.».t < 4'- g ,n,4. d fc'.n.4.4f e i4. 1°) 
k,k f ,q,n 

(B5) 

and find 

U ' 

^U = J^~Y T ( k ' k ''l) e k+q,t e k^ e k'-q,l e k',i (B6) 
k,k' ,q 

where the primed sum is restricted to values of k, k' , k + 
q,k' — q to the domain of the a-type edge state and 

r q) = k +q k k q k B? 

1 - u h ,u k u h ,_u k , 



Appendix C: Background density at an edge 

For convenience, we neglect the spin-degree of freedom 
in the following discussion. Since we do not use explicitly 
spin-dependent Hamiltonians, like the spin-orbit interac- 
tion, the electron spin only leads to a factor of 2. The 
strategy of the following argumentation is to show that 
a half-filled ground state of an aa-ribbon has an elec- 
tron density of exactly | (or 1 with spin) at each site, 
including the sites at the edges. Then, it is shown that 
at half-filling, half of the edge states are occupied. Sub- 
tracting the part of the electron density, originating from 
those filled edge states, from the uniform density leads 
then to the electron density of a state in which all valence 
(conduction) band states are filled (empty) and all edge 
states are empty. 

We consider an aa-ribbon with N y unit cells in the 
transverse direction, described by the Hamiltonian 

N y -i Ny 
H =Y d l,k,A dn ^ B +Y Ukd n,k,A d n-l^B + h.C. 
ri—1 n—1 

. (C1) 

and an odd number N x of unit cells in x-direction, in 
order to exclude k = n. Because of this exclusion, there 
are no states with exactly zero energy. Only eigencnergies 
which are exponentially small in N y exist. The particle- 
hole transformation U = W acts onto the ci-operators 

d n ,k,A — > d n> k,At d n: k,B — > —d n ,k,B- (C2) 

Thus, UHU^ = —H. The density operators dt k s d n ,k,s 
are invariant under U. We now fill exactly half of the 
states 

Ixo) =n©(-e«)4|0> (C3) 

n 

where a^ n creates the eigenstate n with eigenvalue e„ . The 
dual state is |xo) = U \xo) — Y\ n 0( £ n)a} l |0). Because of 
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the particle-hole antisymmetry of the Hamiltonian and 
particle-hole symmetry of the density operators we have 

(Xol d\ k s d n ^ s |xo) = (xo| dl Xs d n , k , s \xo) = \- (C4) 

The states of which |xo) is composed are all valence band 
states plus only the antisymmetric (or only the symmet- 
ric) combinations of the oppositing edge states. This is 
because U turns antisymmetric edge state combinations 
into symmetric combinations. Thus, the part of the elec- 
tron density, \ per site, deriving from the edge state wave 
functions ■i/'o.fc^i s ) is 

where po(n,s) has been defined in Eq. ([7]). It follows 
that the background density is 



PB{n,s) = i - -po(n,s) 



(C6) 



relative to zero filling, or 



p B (n,s) = p B (n,s) - - = --p (n,s) (C7) 



relative to half-filling. 



Appendix D: Comparison to a numerical calculation 

In the main part of this paper, we treated the interac- 
tion as a perturbation to the hopping Hamiltonian Ho. 
Now, we perform numerical lattice-based mean-field cal- 
culations for finite size ribbons. Our numerical calcula- 
tion for an era-ribbon with N unit cells in the transverse 
direction is based on the Hamiltonian 

H = t X] d li,k,A,a( d n,k,B,a + U k d n -l,k,B ,a) + h.C. 
~t e ^2 [dl,k,B : a d O,k,B,a + d* N k A (J d N ^ k . A ^ 



k.a 



u 



\Pn,—q,s,i ) + (pn,q,s,t) Pn-q,s,\\ 



n, q , s 



(Dl) 



with the spin-dependent densities 

/dk 
7^|V , m,fc,(r("-: s )| 2 @ [ e F - e m tT (k)] , 
2-7T 

(D2) 

where ip m ,k{n, s, a) = (0| d nt k, s ,<r |n, k, a) is the wave 
function of the eigenstate to the energy e mjCr (fc). t = 
— 3eV is the hopping amplitude for nearest neighbor hop- 
ping of graphenes 7r-band, t e is the strength of the edge 
gate by which we model a certain class of edge/interface 
properties (see Sect. II), and U is the strength of the 
Hubbard interaction. 



We choose a discrete set of about 4000 k points be- 
tween and 27r in order to approximate the integral in 
Eq. (|D2[) . The Fermi energy is chosen such that exactly 
half of all eigenstates are filled. We start with a density 
that has a small positive magnetization on the one edge 
and a small negative magnetization at the other edge, 
i-e. (p n ,o lS ,a) init = \ + f[j (S n ,oS s ,B - &n,N& s> A), calculate 
the eigenvalue and eigenmodes of (jDll) and from it a new 
set of spin-dependent densities by Eq. (|D2|) . This pro- 
cedure is then iterated until the densities do not change 
anymore, i.e. self-consistence is reached. 
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FIG. 9: (Color online) Comparison of the edge magnetizations 
extracted from the result of the numerical calculation and 
from the effective model. The solid black line is the edge 
magnetization from the effective model. The solid gray lines 
show M" um of an aa-ribbon with different widths (50, 100, 
and 200 unit cells in the transverse direction). The lower 
abscissas of both parts of the figure shows t e /U where the 
band width parameter has been fixed to t e = 0.02eV in Part 
(a) and t e = 0.2eV in Part (b). The corresponding Hubbard 
interaction strength U is shown in the upper abscissas. The 
dashed lines show the polarization of the edge states M° um,es . 

As a quantifier of the solution we choose the edge mag- 
netization. We calculate this quantity for one edge from 
the self-consistent result of the numerical calculation, i.e. 



Ml 



° (Pn 



=0,9=0,5,. 



(D3) 



Note that the M" um is not equal to M e , as defined in Eq. 
([r?j|) . The edge magnetization M e of the effective model 
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only respects the spin-polarization of the edge states. 
M" um , on the other hand, also takes into account the 
polarization of the bulk states near the edges. This ad- 
ditional spin-polarization from the bulk states is small 
for sufficiently weak U (see Fig. Uta)), while for large 
U, the edge magnetization is considerably enhanced (see 
Fig. |U(b)). However, this does not mean that the effec- 
tive edge state model would be insufficient for describing 
the edge magnetism. It only means that the polarization 
of the edge states induces an additional polarization in 
the bulk states via the strong Hubbard interaction which 
increases M" um . In this regime, the edge states are al- 
ready completely polarized. Thus, the difference between 
the effective model and the numerical calculation is only 
quantitative, as desired. 

Further confidence in the validity of the effective model 
can be gained by directly calculating the polarization of 
the edge states, i.e. 

M n„ = [° g|^ , fe , CT (0,5)| 2 e [e F - e . a (fc)] , 

a ~ 

(D4) 

where ipo,k,tr{n, s) is the wave function of the edge state 
with spin a in unit cell n and sublattice site s, with energy 

Af° um,cs is plotted, for different ribbon widths and t e , 
as dashed lines in Fig. [9l Compared to M" um , it is 



much closer to the result of the effective model, as ex- 
pected. Nevertheless, the saturation polarization of the 
edge states is reached already for smaller U. This can be 
interpreted as a back action of the induced polarization 
of the bulk states on the edge states which acts like an 
additional magnetic field. Note that this interpretation 
is in consistence with Ref. 0, where it is found that inte- 
grating out the bulk states in graphene structures leads 
to enhanced effective interaction parameters for the edge 
states. As expected, we observe that such mechanisms 
are only important for large U (see Fig. EJb)) while they 
are absent for small U as can be seen from Part (a) of 
Fig. 1 

Thus, the following physical picture emerges from the 
comparison between the numerics and the analytical 
model: the edge states are primarily responsible for the 
edge magnetization. The polarization of the edge states 
then induces an additional spin-polarization in the bulk 
states if the Hubbard interaction is large enough. This 
additional bulk state polarization further increases the 
edge state polarization so that the saturation is reached 
already for smaller U. The essential approximation, we 
have made in the effective edge state model, is that we 
neglected the enhancement of the effective interaction by 
the bulk states. This enhancement can be easily reintro- 
duced into the model, if desired. 



1 K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 
Zhang, S. V. Dubonos, I. V. Grigorieva, and A. A. Firsov, 
Science 306, 666 (2004). 

2 A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 

3 M. Y. Han, B. Ozyilmaz, Y. Zhang, and P. Kim, Phys. 
Rev. Lett. 98, 206805 (2007). 

4 F. Sols, F. Guinea, and A. H. Castro Neto, Phys. Rev. 
Lett. 99, 166803 (2007). 

5 M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, 
J. Phys. Soc. Jpn. 65, 1920 (1996). 

6 J. Jung and A. H. MacDonald, Phys. Rev. B 79, 235433 
(2009). 

7 K.-I. Sasaki and R. Saito, J. Phys. Soc. Jpn. 77, 054703 
(2008). 

8 Y.-W. Son, M. L. Cohen, and S. G. Louie, Phys. Rev. Lett. 
97, 216803 (2006). 

9 L. Pisani, J. A. Chan, B. Montanari, and N. M. Harrison, 
Phys. Rev. B. 75, 064418 (2007). 

10 H. Feldner, Z. Y. Meng, A. Honecker, D. Cabra, S. Wessel, 
and F. F. Assaad, Phys. Rev. B 81, 115416 (2010). 

11 S. Dutta, S. Lakshmi, S. K. Pati, Phys. Rev. B 77, 073412 
(2008). 

12 T. Hikihara, X. Hu, H.-H. Lin, C.-Y. Mou, Phys. Rev. B 
68, 035432 (2003). 

13 B. Wunsch, T. Stauber, F. Sols, and F. Guinea, Phys. Rev. 
Lett. 101, 036803 (2008). 

14 E. H. Lieb, Phys. Rev. Lett. 62, 1201 (1989). 



E. Lieb and D. Mattis, Phys. Rev. 125, 164 (1962). 

M. J. Schmidt and D. Loss, arXiv 2009. 

K. Sasaki, S. Murakami, and R. Saito, Appl. Phys. Lett. 

88, 113110 (2006). 

T. Giamarchi, Quantum Physics in One Dimension, Ox- 
ford University Press, 2004. 

C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

The movies can be downloaded from 

|http: / /quantumtheory.unibas.ch/schmidtm/edge jnagnetism.html 
L. Bartosch, M. Kollar, and P. Kopietz, Phys. Rev. B 67, 
092403 (2003). 

Joaquin E. Drut and Timo A. Lahde, Phys. Rev. B 79, 
165425 (2009). 

M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, 
J. Phys. Soc. Jpn. 65, 1920 (1996). 

M. Kinza, J. Ortloff, and C. Honerkamp arXiv:1002.2513vl 
(2010). 

Because the support of the a edge state in the Brillouin 
zone of a GNR is k £ the maximum momen- 

tum difference at the Fermi levels fci?jj — hp,L < ^g-, so 
that first order umklapp processes are not allowed. For /3 
edge states the situation is different and for some filling 
fractions, umklapp scattering may well play an important 
role. 

We exclude interactions that are strong enough to create a 
bulk gap in graphene 2 ^ thus breaking sublattice symmetry 
spontaneously. This may be done by placing the structure 
on a Si02 substrate. 



We use the notation from Ref. [Isl . 



